Nonlinear equation for anomalous diffusion: unified power-law and stretched 

exponential exact solution 
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The nonlinear diffusion equation ^ — j^^p lo ci,iicii_y iic±c, wii^io ^ — 
d, 6 and v are real parameters. This equation unifies the anomalous diffusion equation on fractals 
{u — 1) and the spherical anomalous diffusion for porous media (6 = 0). Exact point-source solution 
is obtained, enabling us to describe a large class of subdiffusion {6 > (1 — i')d), normal diffusion 
{6 — {1 — u)d) and superdiffusion {9 < {1 — v)d). Furthermore, a thermostatistical basis for this 
solution is given from the maximum entropic principle applied to the Tsallis entropy. 
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One of the most ubiquitous processes in nature is the 
diffusive one. In this context, anomalous diffusion has 
awakened great interest nowadays, in particular in a va- 
riety of physical applications. A representative set of 
such applications of current interest is surface growth 
and transport of fluid in porous media [Q, diffusion in 
plasmas ||] , diffusion on fractals Q , subrecoil laser cool- 
ing jj], CTAB micelles dissolved in salted water [||, two 
dimensional rotating flow |^] and anomalous diffusion at 
liquid surfaces [Q. The anomalous diffusive process is 
commonly characterized from the mean-square displace- 
ment time dependence, (r^) oc , with cr ^ 1, i.e., we 
have superdiffusion for a > 1 and subdiffusion for cr < 1. 
For a system which presents anomalous spreading, it is 
generally associated to a non-gaussian space-time distri- 
bution, like power-law or stretched exponential. In this 
framework, it is desirable to incorporate, in a unifled way, 
these two behaviours, since it enables us to describe a 
wide class of diffusive processes. The present work is 
dedicated to giving such unified description. 

Power-law or stretched exponential distributions arise 
naturally from generalizations of the d-dimensional dif- 
fusion equation 



(1) 



with p = p{x,t), X = {xi,X2,. ■ ■ ,Xd), A ^ 
X]n=i ldx\, and D being the diffusion coefficient. The 
nonlinear equation 



(2) 



is just one of these generalizations, where is a real pa- 
rameter. Eq. ^ has been employed to model diffusion in 
porous medium (see Ref. [Q and references therein) and 
in connection with generalized Tsallis statistics [p|-p^. 
Another important kind of anomalous diffusion, in a 



tridimensional space, is related to turbulent diffusion in 
the atmosphere and is usually described by [0| 



dp 
dt 



= V • (KVp) , 



(3) 



where K oc r^^^ (r = |a;|). In a more general case, we 
consider Eq. (|3|) in a d-dimensional space with K oc r~^, 
where is a real parameter. Thus, V • {K\7) is propor- 
tional to r"(''"i)^r'*~^~''^ -f A/r^^'^ {A is an operator 
depending on the angular variables), and consequently 
A = r~^'^~^'> ■^r'^^^^^ is the radial part to be consid- 
ered in the study of the spherical symmetrical solutions 
of Eq. (^). In this context, when d is interpreted as frac- 
tal dimension in an embedding iV-dimensional space, the 
equation 



dp 
dt 



DAp 



(4) 



has been used to study diffusion on fractals 
Here we are going to propose the equation 

| = v.(W), 



(5) 



as unification of Eqs. (||) and (^). In fact, Eq. (||) reduces 
to the correlated anomalous diffusion (0) if K = D, and 
to the generalized Richardson equation (0) if = 1. The 
present study is mainly addressed to the point-source so- 
lution of Eq. (^, because it contains, as particular cases, 
a(n) (assimptotic) power-law and a stretched exponen- 
tial. In this away, we focus our attention on the radial 
equation 



(6) 



Using this equation instead of Eq. (||) enables us to an- 
alyze cases with noninteger d, so we can relate d with 



1 



a fractal dimension. Therefore, in the following discus- 
sion, we are going to consider c? as a nonnegative real 
parameter. 

In order to motivate the ansatz to obtain an exact time- 
dependent solution for Eq. , we recall the correspond- 
ing solutions for equations Eqs. (|l|), (||), and (||). The 
time-dependent point-source solution for Eq. p) is 



Po 



■ exp 



(7) 



where the normalization p{r,t)r'^ ^dr = po and 

the d-dimensional solide angle fi^ = 27r'^/^/r((i/2) have 
been used. From Eq. (Q) we can easily obtain the Ein- 
stein formula for the Brownian montion, i.e., (r^) — 
2dDt. 

The analogous solution for Eq. (||) is jS- 10 



pir,t) = [l-{l~q)/3,{t)r^ 



/Zi{t) 



(8) 



where q = 2 - Zi{t) oc t'^/i^+dii-g)] ^nd (3i{t) oc 
^-2/[2+rf(i-(/)] jg important to emphasize the short 
or long tailed shape of Eq. (^) , when compared with the 
normal diffusion (limit q 1). When q < 1 we have 
p{r,t) = for 1 — (1 — q)l3i{t)r'^ < 0, giving the short 
tailed behaviour for p{r,t). On the other hand, when 
(? > 1, the asymptotic power-law behaviour for the so- 
lution 7.-2/(9-1)^ shows that p{r,t) is a long tailed 
function. This short or long tailed behaviour for p(r, t) 
reflects directly on the mean-square displacement, lead- 
ing to (r^) oc t2/[2+(i(i-g)]_ Again compared with the 
usual diffusion, q = 1, we have a superdiffusion (subdif- 
fusion) for q > l(g < 1). 

The fundamental solution of Eq. (|^) is the stretched 
exponential 



p(r,t)=exp(-/32(t)/+2)/Z2(i), 



(9) 



which Z2{t) oc and (32{t) oc t~^, presenting a 

short (long) tailed behavior for > (0 < 0). Fur- 
thermore, the mean-square displacement behaviour is 
(r2) oc t2/(9+2)_ Thus, for 6* > (6* < 0) we have subdif- 
fusive (superdifussivc) regime. 

Note that Eqs. (Q), (g|) and (H) can be interpolated 
if we employ a generalized stretched Gaussian function. 



G 



{x)^[l-{l-q)\x\^f^'''^ or G(,^,)(x) 



when 1 — (1 — < 0, with q and A being real 

parameters. In this direction, our ansatz to solve Eq. (|^) 
is 



p(r,t)= [l-{l-q)(3{t)r 



Al — 



(10) 



or p(r, i) ^ if 1 - {l-q)l3{t)r^ < 0,where P{t) and Z{t) 
are functions to be determined. By using this ansatz in 
Eq. (§) we verify that l3{t) and Z{t), with A = 61 + 2 and 
q — 2 ~ ly, obey the equations 



dZjt) 

dt 
dpjt) 

dt 



= DX{2-q)d/3{t)Zi{t) 

= -DX\2~q)p^{t)Z^-\t). 



(11) 



The solutions of these nonlinear differential equations, 
that lead Eqs. (0), (||) and (||) as limit cases of Eq. (|l^), 
are 



l3{t) = A t ^+"(1-9) 
where 



Z{t)^Bt' 



(12) 



^= {7'-^ [D\{2-q){X + d{l-q))]] 



B 



which 



DX{2-q){X + dil-q)) 



(13) 



2^ £11) 

Apor(i) 



(9-i)*r(^) 



r.d-1, 



(14) 



The normalization condition, p{r,t)r'^^^dr — po, 

employed in the above calculation can only be satisfied if 
A > and X+d{l—q) > 0. From these conditions over the 
parameters q and A, we verify that the exponents in f3{t) 
and Z{t) are respectively negative and positive. In addi- 
tion to the normalization condition, the restriction q < 2 
is necessary for p{r, t) to be real. In the following, we 
assume that the parameters obey the above restrictions. 
Of course, by setting the appropriate limits of parame- 
ters and V, or equivalently A and q, the solutions (|^), 
(^ and (p|) are recovered, giving the full expression for 
(3i{t), Zi(t), [32{t) and Z2{t). 

By using the above solution we can calculate the mean 
value of ; it is 

(0 = ft^^^^^^^=C„t^, (15) 



where 



/p p{r,t)r'^ ^dr 

r(^-^) 



Y(d+<i\ (g-i)Tr(^-f) 



r(f) 



r(4+T^+i) 



(i-9)-r(^+T^+i) 



('Z>1) 



('Z<1) 



with A given by Eq. (|T^). When g < 1, the mean value 
(r") always exist. On the other hand, the existence of 
(r") for q > 1 impose a further restriction over the pa- 
rameters: A + d{\ — q) > a{q — 1). 

To decide if the diffusion is anomalous or normal, 
we consider ([l5| ) with a — 2. In this way, we have 
< >oc r with a = 2/[2 + 9 + d{v - 1)]. Thus, the 



2 



condition for normal diffusion, cr = 1, can be satisfied 
even when p does not obey Eq. (0), i.e., 9 = d(l — v) 
with and v ^ \. In this case, we can also verify that 
the anomalous diffusive regime induced by 7^ is com- 
pensated by a convenient one with v ^ 1. Furthermore, 
this competition between and v values can lead to a 
subdiffusion (cr < 1) if 6* > — i^) or a supperdiffusion 
{a > 1) \i 9 < d{l — v). This classification is illustrated 
in Fig. (I) for d = 3. 

In the following, we discuss the consequences of the 
above classification on the p shape. From the solution 
(^, for porous medium, and the solution (|^), for diffu- 
sion on fractals, we can see that the supcrdiffusive (subd- 
iffusive) regime is associated to long (short) tail of p(r, t) 
when compared with Gaussian (0). However, this conec- 
tion is not valid in general. To illustrate the relation 
between regime of diffusion and tail behaviour of p(r, t), 
we consider Fig. (|^). In this figure, we plot p{r,t) given 
by Eq. (p^ ) versus r to some values of 9 and v, subject 
to the restriction 9 = d{l — (the normal diffusion line 
indicated in Fig. (|l])). In this case we observed short and 
long tail behaviours compared with the Gaussian one. 



where Zq = /^[l - (1 - q)P\x\^]^/'^^-iUx and /? is re- 
lated to the Lagrange multipliers. 




FIG. 1. The diagram indicates the diflusive regime related 
to Eq. ^ in terms of its parameters 6 and v to d — 3. 
Thus, by using (r^) oc f2/[e+2+d{..-i)] ^^^^ y^^^^ 

classified the subdiffusive {0 > {l — i')d), normal (6' — (1— 
and superdiffusive {9 < (1 — i')d) regimes. The forbidden 
region refers to the region of parameters where (r^) does not 
exist (diverges). 

To conclude our discussion about Eq. (||) and its ra- 
dial time-dependent solution (|o|) , we present an entropic 
basis for this solution. This basis is motivated by the 
Tsallis generalized statistical mechanics p"^-p^, where 
the TsaUis entropy |l| 5, = (1 - T,Y=iPj)/i^ " <?) plays 
a central role with {pj} being the probabilities for the 
W states of the system, and q ^ TZ being the Tsallis 
index (by taking the limit g — > 1 we recover the usual en- 
tropy iS*! = — J2j=i P] I'^Pj)- To understand the entropic 
basis, for simplicity, let us consider the maximization of 
Sq = {l — J'^^p{x)'^dx)/{q — l) subject to the constraints 

@S iZo = 1 and \x\^p{x) dx = Uq. This 

maximization leads to p(2:) = [l — {l — q)(3\x\^]^^^-^~'^''/Zq, 




FIG. 2. By considering the normal diffusion, 6 — {1 — v)d, 
the shape of p(r, t) with d = 3 and t = 5 is illustrated in three 
cases: short t&i\{6 = —0.6 and v = 1.2), Gaussian [9 = and 
V = 1), and long tail [9 = 0.6 and v = 0.8). Inset plot: detail 
of tail behaviour in the three cases above. 
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